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摘 要 : 由 于 采 后 处 理 过程 中 脐橙 保鲜 剂 抑 霉 唑 易 通过 果皮 渗 进 果肉 中 残留 ， 不 慎 食 用 后 会 对 人 体 产生 危 


害 。 因 此 ， 本 研究 探索 一 种 基于 表面 增强 拉 曼 光谱 技术 (Surface-Enhanced Raman Spectroscopy, SERS) 的 


脐橙 果皮 中 抑 霉 唑 残留 的 快速 检测 方法 。 首 先 对 SERS 检 测 条 件 进行 优化 ， 分 别 确定 了 最 优 的 检测 条 件 为 反 


应 时 间 2 min， 金 胶 加 入 量 400 uL, NaBr EX H 


有 解 质 溶液 量 加 入 量 为 25 pL。 基 于 以 上 最 优 检测 条 件 ， 以 自 
JH TE CART tie) = FEI (Adaptive Iterative Reweighted Penalized Least Squares, air PLS), air PLS+ 归 一 化 、 
air PLS+ 基 线 校正 、air PLS+ 一 阶 导 数 、air PLS+ 标 准 正 态 变量 (Standard Normal Distribution, SNV) 和 air 


PLS+ 多 元 散射 校正 (Multiplicative Scatter Correction，MSC) 处 理 后 的 6 组 光谱 数据 为 研究 对 象 ， 分 别 采 用 这 
6 种 光谱 预 处 理 法 建立 支持 向 量 回归 (Support Vector Regression, SVR) 模型 并 对 预测 性 能 进行 比较 后 发 现 ， 
air PLS 方 法 所 建立 模型 的 预测 集 相关 系数 (Coefficient of the Determinant for the Prediction Set, R,) 最 大 ， 预 
测 集 均 方 根 误差 (Root-Mean-Square Error of Prediction, RMSEP) 最 小 。 对 光谱 数据 进行 主 成 分 分 析 (Prin- 
cipal Component Analysis, PCA) 特征 提取 ， 选 择 前 7 个 主 成 分 得 分 作为 SVR 预测 模型 的 输入 值 。 采 用 SVR、 
多 元 线性 回归 (Multiple Linear Regression, MLR) 和 偏 最 小 二 乘 回归 (Partial Least Squares Regression, 
PLSR) 三 种 建 模 方法 分 析 比 较 其 对 应 的 预测 性 能 ， 其 中 SVR 模 型 的 预测 集 , 可 高 达 0.9156， 预 测 集 RMSEP 
为 4.8407 mg/kg， 相 对 标准 偏差 (Relative Standard Deviation, RPD) 为 2.3103， 表 明基 于 SVR 算法 对 脐橙 表 
面 抑 霉 唑 残留 的 预测 值 越 接 近 实 测 值 ， 越 能 有 效 提高 模型 预测 准确 性 。 试 验 结果 表明 ， 利 用 SERS 结 合 PCA 


及 SVR 建 模 ， 可 实现 对 脐橙 果皮 中 抑 才 唑 残留 的 快速 检测 。 
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动 了 当地 经 济 的 发 展 ， 深 受 当地 政府 重视 ,该 产 
地 也 成 为 了 全 国 最 大 的 脐橙 生产 区 O. PEE 
藏 过 程 中 易 受 病原 菌 侵害 ,导致 其 腐烂 ， 造 成 经 
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济 损失 2” 。 目 前 ， 采 用 化 学 保鲜 剂 是 柑橘 类 水 果 
储存 中 病害 防治 相对 成 熟 、 效 果 较 好 的 保鲜 方 
st. MÆ (Imazalil, IMZ) 和 咪 鲜 胶 等 化 学 保 
鲜 剂 主要 用 来 防治 脐橙 的 青 霉 病 和 绿 霉 病 O, 
但 该 类 保鲜 剂 长 期 使 用 容易 产生 耐 药 性 ， 进 而 引 
起 保鲜 效率 的 折 损 。 不 少 商 贩 为 了 达到 保鲜 效 
果 ， 在 贮藏 过 程 中 加 大 保鲜 剂 的 使 用 剂量 ， 造 成 
脐 检 中 抑 考 唑 过 量 残 留 。 而 且 ， 即 便 是 少量 的 抑 
霉 唑 残留 也 会 使 动物 和 人 类 的 肝脏 损伤 ， 甚 至 可 
能 致癌 “。 为 此 ， 中 国 规定 了 脐橙 中 抑 霉 唑 保 
鲜 剂 残留 的 最 大 限量 为 5 mg/kg, OW T HIRI 
橙 的 食用 安全 品质 ， 检 测 脐橙 果皮 中 抑 霉 唑 的 残 
留 量 具有 重要 的 实际 意义 。 

目前 常用 的 抑 考 唑 残留 检测 方法 主要 有 人 气相 
色谱 法 (Gas Chromatography, GC) ^", Xi 
相 色 谱 法 (High Performance Liquid Chromatogra- 
phy, HPLC) ^", iti? (Mass Spectrometry, 
MS)"” 和 超 高 效 液 相 色谱 法 (Ultra-Perfor- 
mance Liquid Chromatography, UPLC) |" 等 。 以 
上 方法 主要 优点 在 于 再 现 性 好 、 精 确 性 和 敏感 性 
高 ， 但 所 需要 的 仪器 设备 较 昂 贵 ”"， 试 验 操作 
繁琐 ， 耗 时 长 ， 无 法 实现 对 抑 霉 唑 残留 量 的 快速 
检测 ， 在 实际 应 用 中 受到 一 定 的 限制 。 因 此 ， 有 
必要 探索 一 种 可 以 快速 检测 脐橙 果皮 中 抑 考 唑 残 
留 的 方法 。 

表面 增强 拉 曼 光谱 技术 (Surface-Enhanced 
Raman Spectroscopy, SERS) 是 一 种 当 入 射 激光 
照射 到 金 或 银 等 粗糙 金属 表面 时 ， 在 其 表面 会 激 
发 出 等 离子 体 并 在 粗糙 金属 粒子 的 覆盖 区 域内 产 
生 共 振 ， 使 其 表面 产生 比 普 通 拉 曼 信号 强 4~11 
个 数量 级 的 信号 强度 的 技术 ， 该 方法 具有 方便 快 
捷 、 灵 敏 度 高 、 稳 定性 好 、 操 作 简 单 和 可 适用 于 
现场 分 析 的 优点 o Zhao $ U R HH SERS f 
术 制 备 了 柔性 衬 底 ， 对 苹果 表面 抑 霉 只 残留 进行 
原 位 检测 ， 在 苹果 上 检测 到 的 最 低 浓度 可 达 
0.073 mg/kg。 该 研究 表明 ，SERS 方法 具有 较 高 
的 灵敏 度 。 抑 霉 唑 在 水 果 采 后 保鲜 的 过 程 起 到 非 
Wo BEM ED. Ayah Se h PBT EAR SZ I sO GF 


问题 对 人 体 造 成 伤害 ， 采 用 SERS BOR FS DU or RE 
果皮 中 抑 霉 唑 残留 具有 重要 意义 。 本 人 研究 主要 以 
脐橙 为 对 象 ， 采 用 自 适 应 迭代 惩罚 最 小 二 乘法 
(Adaptive Iterative Reweighted Penalized Least 
Squares, air PLS) 对 含 不 同 浓度 的 抑 霉 唑 的 脐橙 
果皮 的 SERS 光谱 进行 预 处 理 ， 并 选择 前 7 个 主 
成 分 得 分 值 作为 支持 向 量 回归 (Support Vector 
Regression, SVR) 模型 的 输入 值 ， 建 立 SVR 预 
测 模型 对 脐橙 果皮 中 的 抑 霍 唑 残留 量 进行 预测 。 


2 材料 与 方法 


2.1 仪器 与 设备 

便携 式 拉 曼 光 谱 仪 系统 主要 包括 QE65Pro 型 
拉 曼 光谱 仪 (美国 海洋 光学 公司 )、785 nm 激光 
器 、 光 纤 和 采样 附件 等 ,DXR™ 显 微 拉 曼 光谱 仪 
(美国 赛 默 飞 世 尔 科 技 有 限 公 司 ); 磁力 搅拌 器 
(LC-DMS-H， 上 海 力 展 邦 西 仪器 科技 有 限 公 
HD; 数控 超声 波 清 洗 器 (KQ-500DE 型 ， 昆 山 
市 超声 仪器 有 限 公司 ); 电子 天 平 (FA2004， 上 
海 玛 宇 恒 平 科学 仪器 有 限 公 司 ); 实验 室 超 纯 水 
机 (美国 TeLedyne 水 质 公 司 )。 


2.2 材料 与 试剂 


本 研究 所 用 脐橙 购 于 江西 农业 大 学 水 果 店 。 
研究 所 用 试剂 包括 三 水 合 四 氧 金 酸 (HAuCl 
3HO>49.0%， 西 格 玛 奥 德里 奇 贸 易 有 限 公司 ) , 
抑 霉 只 乳油 (500 g/L, I ETEA RA 
HJ), MEERE m (IMZ， 纯 度 >99%， 阿 拉 
TO. ， 柠 檬 酸 三 钠 (CH,NaO,2H,O, 99%, P 
陇 化 工 股 份 有 限 公 司 ) ， 澳 化 钙 (CaBr,, 298%, 
阿拉 丁 )， 省 化 钠 (NaBr, 99%, BHT), AE 
钙 (CaClL， 分 析 纯 ， 纯 度 98%) ， 和 氧化 钠 (Na- 
Cl1，99.5%， 西 陇 化 工 股份 有 限 公 司 ) ， 以 及 实验 
室 自制 超 纯 水 。 


2.3 样本 制备 


(1) 抑 霉 唑 原液 配制 。 取 100 uL AY Hl He eE 
乳油 (500 g/L) 置 于 100 mL ABH, MA 
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62.5 mL 的 超 纯 水 均匀 混合 ， 得 到 800 mg/L B ffl 
霉 唑 原液 。 

(2) 抑 霉 唑 工作 溶液 制备 。 取 不 同体 积 的 乳 
il (500 g/L) 用 超 纯 水 稀释 成 不 同 质量 浓度 的 
METERA (250—500 mg/L). 

(3) 样本 制备 。 将 脐橙 的 果皮 用 超 纯 水 洗 
净 ， 蚊 干 后 将 果皮 切 成 若干 个 面积 为 1.5 cemX 
1.5 cm WAREZ) ER, RET Pe FE (a, FA 
位 : g)o 

在 条 件 优 化 试验 中 ， 将 切 好 的 脐橙 小 块 浸 
泡 于 抑 霉 哗 深 液 (800 mg/L) 中 ，90 s 后 捞 出 
BUT. FOR (b， 单 位 : g)， 得 到 含 
40.34 mg/kg 抑 霉 唑 的 脐橙 果皮 样本 。 

在 定量 试验 中 ,将 切 好 的 65 个 脐橙 小 块 浸 
FMM (250—500 mg/L) 中 ，90 s 后 捞 
出 景 干 ， 再 对 其 称 重 (b, Hu: g)， 得 到 含 抑 
霉 唑 的 脐橙 果皮 样本 (1.43~40.81 mg/kg). A 
浸泡 在 抑 考 唑 中 的 脐橙 果皮 样本 为 空白 脐 覆 果 皮 
样本 。 其 中 ， 在 条 件 优化 试验 部 分 与 定量 试验 部 
分 中 ，a 为 脐橙 果皮 的 初始 重量 ，b 为 浸泡 抑 霉 
只 溶液 且 晾 干 后 的 重量 ; b-a 可 得 到 脐橙 果皮 中 
MERKET, g; 设 c 为 抑 霉 唑 洲 液 的 浓度 ， 
mg/kg, J| (b-a) X c J PrE R Rz THU A SIC s c ifc 
{AEM AY Hat, mg; (b-a) X c/b W WEAR ESR 
皮 上 实际 吸收 的 抑 霉 唑 的 浓度 ，mg/kg。 

(4) 金 胶 制 备 。 取 3 mL 浓度 为 1% 的 
HAuCL, 溶 液 于 150 mL 的 烧杯 中 ， 再 慢 慢 加 入 
47 mL 超 纯 水 ,混合 均匀 ， 将 混合 液 放 至 磁力 搅 
拌 器 中 加 热 到 沸腾 。 沸 腾 后 立即 加 入 2 mL 浓度 
为 1% 的 柠檬 酸 三 钠 溶液 ， 并 搅拌 8 min Ja, 8 
却 至 室温 形成 金 胶 溶 液 备 用 。 

(5) 电解 质 溶液 制备 。 取 0.1 g CaBr,、 0.55g 
CaCl,、5.14 g NaBr 和 2.9 g NaCl 分 别 置 于 棕色 的 
50 mL 容量 瓶 中 ， 先 加 入 少量 水 溶解 ， 再 定 容 至 
50 mL， 可 分 别 得 到 0.01 mol/L 的 CaBr, 溶液 、 
0.1 mol/L CaCL7& X. 1 mol/L NaBr 和 1 mol/L 
NaCl 溶 液 。 


2.4 拉 曼 光谱 采集 


抑 霉 只 标准 品 拉 曼 光谱 采集 。 取 适量 的 
抑 霉 唑 标准 品 于 载 玻 片 上 ， 用 DXRY 显 微 拉 
曼 光 谱 仪 采集 其 拉 曼 光谱 。 本 研究 主要 分 析 
400~1800 cm ' 波 段 的 拉 曼 光谱 。 

定性 试验 中 脐橙 果皮 SERS 光谱 采集 。 将 
400 wL 纳米 金 胶 和 25 uL NaBr 溶 液 充分 混合 后 喷 
到 含 40.34 mg/kg MER AY RE NER D, FER 
SiMe) FIRM 2 min 后 ， 采 集 其 SERS 光谱 。 
每 个 优化 条 件 下 ， 设 置 5 个 平行 样 ， 取 其 平均 值 
为 该 条 件 下 的 拉 曼 光谱 。 

定量 试验 中 脐橙 果皮 SERS 光谱 采集 。 将 
400 uL 纳米 金 腕 和 25 uL NaBr 溶 液 充分 混合 后 喷 
Bl T d AE ME AY Hb RE 7 ER (1.43—40.81 mg/kg) 
E, REPE SERS. BETSY BEF BD 
样本 采集 5 条 拉 曼 光谱 作为 平行 样 ， 取 其 平均 值 
作为 该 样本 的 拉 曼 光谱 。 

其 中 ，DXR™ 显 微 拉 曼 光 谱 仪 的 参数 设置 为 : 
激光 能 量 5mW， 激 发 波长 780 nm， 物 镜 的 放大 倍 
数 选 择 10 倍 ， 采 集 曝 光 时 间 10 s， 预 览 采 集 时 间 
10 s， 样 品 曝光 10 次 ， 背景 曝光 116 次 。 便 携 式 
拉 曼 光谱 仪 的 参数 设置 为 : 积分 时 间 为 60s， 平 
均 次 数 为 2， 平滑 度 为 1。 


2.5 SERS 检 测试 验方 法 


2.5.1 定性 检测 

3B iE EE Be as A IB RE AR Bey MI UP] JURE 
果皮 (40.34 mg/kg) AY SERS 光谱 ， 得 到 983 
All 1047 cm" JE BS VEN Fe E WS ol] JE ESR A E 
否 残 和 留 抑 霉 哗 。 因 此 ， 本 研究 通过 比较 983 和 
1047 cm 波段 处 的 SERS 信号 强度 来 确定 最 佳 的 
检测 条 件 。 

为 研究 不 同 反应 时 间 对 SERS 强度 的 影响 ， 
将 400 uL 4x fé 25 kL 的 NaBr 溶 液 混合 均匀 后 
喷洒 到 含 40.34 mg/kg ME BS RE RE AS b. OP 
别 采集 反应 时 间 为 0、2、4、6 和 8 min 时 的 
SERS 光谱 。 通 过 比较 983 和 1047 em" JE Et Ah AY 
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SERS fri v8 HE 53 oz JW EST [EU AY IG AOR fü iat HD n 
分 子 与 金 胶 的 最 佳 反应 时 间 。 

为 研究 不 同 金 胶 加 入 量 对 SERS 强度 的 影 
响 ， 分 别 将 80、240、400、560 和 720 uL 的 金 


曼 光 谱 频 率 ， 在 计算 时 采用 DFT 方法 中 的 
B3LYP/6-31G (d, p) 基 组 。C、H、CI、N 和 O 
原子 也 采用 6-31G (d, p) 基 组 ， 在 Gaussian 09 
中 显示 计算 结果 。 对 拉 曼 位 移 用 该 机 组 相对 应 的 


胶 与 25 uL 的 NaBr 溶 液 混合 均匀 后 喷 酒 到 合 
40.34 mg/kg TIE WE AY TREE AS E. 52 min 后 
采集 其 SERS 光谱 。 通 过 比较 983 和 1047 cm X 
段 处 的 SERS 信号 强度 与 金 胶 加 入 量 的 关系 来 确 
定 最 佳 的 金 胶 量 。 

为 研究 不 同 电 解 质 溶液 的 加 入 量 对 SERS 强 
度 的 影响 ， 将 400 pL 金 胶 分 别 与 25 uL 的 CaBr, 
WR. Cac RW, NaBr WA NaC RORE 
均匀 后 喷洒 到 含 40.34 mg/kg 抑 霉 唑 的 脐橙 样本 
上 ， 反 应 2 min 后 采集 其 SERS 光谱。 通过 比较 
983 和 1047 cm Wk Be Lb AY) SERS 信和 号 强度 与 电解 
质 溶液 (CaBr,, CaCl,, NaBr 和 NaCl) 的 直方 
图 来 确定 最 佳 的 电解 质 溶液 。 


校正 系数 (0.961) 校正 。 其 中 基 组 为 该 电子 体 
系 内 轨道 的 数学 描述 , “6” 解 释 为 内 层 轨 道 的 函 
数 数量 ; 而 一 个 基 组 中 含 3 个 高 斯 函数 (Gauss- 
ian Type Orbital, GTO) ， 每 一 个 Slater PK 2 
(Slater Type Orbital, STO) 由 6 个 GTO 展开 ， 
“31” 解 释 为 两 组 STO 描述 的 价 层 的 轨道 ;“d” 
和 “p” 分 别 是 对 分 子 中 的 C、Cl、N、O 原 子 和 
HH 原子 做 的 极 化 函数 I 


3 结果 与 分 析 


3.1 HEERE DHS KAM 
抑 霉 唑 的 分 子 式 为 CH,CLNO， 由 1 个 烯 


为 了 研究 不 同 NaBr 加 入 量 对 SERS 强度 
的 有 影响， 将 400 pL Az Hé) O 10, 15, 20, 
25. 30 uL 的 NaBr 深 液 混 合 均匀 后 喷洒 到 含 
40.34 mg/kg MERERI GREE AS E. 5 2 min 后 
采集 其 SERS 光谱 。 通 过 比较 983 和 1047 om" YK 
段 处 的 SERS 信号 强度 与 NaBr 加 入 量 的 折线 图 来 
确定 最 佳 的 电解 质 涂 液 。 
2.5.2 ”定量 检测 

将 400 uL 4: Iz 5j 25 uL AY NaBr IMIG 
35] ^] Jes Wit WG 80 xr Ug ME B Hop RETE (1.43~ 
40.81 mg/kg) E, 反应 2 min 后 采集 其 拉 曼 光谱 。 
得 到 65 个 含 抑 霉 唑 的 脐橙 样本 ， 从 中 随机 挑选 
35 个 建立 预测 模型 ， 剩 余 的 30 个 脐橙 样本 则 用 
作 预 测 集 ， 用 来 评估 该 模型 的 预测 性 能 。 


2.6 MBM ths tists 


为 对 比 理论 数据 结果 和 试验 结果 之 间 的 一 臻 
性 ， 采 用 密度 泛 函 理论 (Density Functional The- 
ory, DFT) 计算 方法 对 抑 考 唑 分 子 的 振动 归属 
进行 理论 运算 分 析 。 

Tp SEMPER 3D 4) T 258192 H] Gauss View 5.0 构 
itt, Gaussian 09W $k (FITS BRE Bek 2) T- B 4v. 


Vj ASE 1 SACRE RE ZA. HLA f 
5| B4] Sg MEY 3-58 gun 1 所 示 ， 主 要 包括 C、 
C-H, C=C-H, C-O, N=C-N, C-N, N-C-H 和 
C-N-C 5& XE H] 


a 


a 
图 1 优化 的 抑 霉 唑 分 子 结构 
Fig. 1 Optimized molecular structure of IMZ 

图 2 (a) AN BYE Hh AY TIS ME hr ETE, 
图 2 (b) SE Pa ill E B) TIRE E p PE a FH] ore 
谱 。 通 过 对 比 两 者 的 特征 峰 归属 可 发 现 ，1176 和 
1378 cm 波段 处 特征 峰 的 拉 曼 信号 强度 较 弱 。 
这 可 能 是 由 于 这 两 个 位 移 处 分 别 在 茶 环 上 有 微弱 
的 面 外 弯曲 和 面 内 弯曲 振动 所 引起 。 结 合 表 1 和 
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图 2 可 得 ， 抑 霉 唑 标准 品 在 690 cm "波段 处 的 特 ” ” 拉 曼 光谱 的 特征 峰 匹 配 可 发 现 ， 最 小 偏差 为 0， 


征 峰 是 由 于 茶 环 弯曲 振动 、CHH 面 内 摇摆 、C-C ”最 大 偏差 为 10。 这 可 能 是 由 于 实际 测量 拉 曼 光谱 


面 外 弯曲 、 二 毛茶 基 面 外 弯曲 共同 作用 产生 的 ; ”的 标准 品 为 固体 粉末 ， 其 之 间 存 在 相互 作用 力 ， 
983 cm "波段 处 的 振动 可 能 是 由 于 茶 环 弯曲 振动 、 ”而 理论 计算 时 的 形态 为 理想 的 气态 单 分 子 ， 忽 略 
二 毛茶 基 上 C=N-C 和 C-N-C 与 C=C-N 的 弯曲 振 ”了 分 子 间作 用 力 。 基 于 此 ， 理论 计 算 的 抑 考 唑 拉 
动 共 同 作 用 产生 ; 1041 em Ab f e e 
振动 、C=C-H 面 外 弯曲 振动 和 CH, 面 内 弯曲 振动 确 性 ， 可 以 将 实测 的 抑 霉 唑 标准 品 拉 曼 光谱 作为 
共同 作用 产生 。 图 2 中 理论 计算 光谱 与 标准 品 的 “ ”后续 分 析 膀 橙 表面 抑 霉 唑 残留 的 参考 依据 。 

表 1 抑 雷 唑 理论 计算 和 标准 品 拉 曼 特征 峰 归 属 


Table 1 Raman peak attribution analysis of IMZ theoretical calculation and standard 


由 C-HR ” 曼 光 谱 很 好 地 印证 了 实测 的 抑 霉 哗 拉 曼 光 谱 的 正 


标准 品 拉 曼 特征 峰 /cm” ”理论 拉 曼 特征 峰 /cm” 主要 特征 峰 归属 
690 686 茶 环 弯曲 振动 .CH, 面 内 摇摆 ,、C-C 面 外 弯曲 二 氧 茶 基 面 外 弯曲 
983 977 茶 环 弯曲 振动 .二 毛茶 基 上 C=N-C 和 C-N-C 与 C=C-N 的 弯曲 振动 
1041 1051 C-H 剪 式 振动 .C=C-H 面 外 弯曲 振动 .CH, 面 内 弯曲 
1176 1177 C-C fll C-O 伸缩 振 动 .C-H 弯 曲 振动 . 葵 环 面 外 弯曲 
1378 1372 N=C-N 剪 式 振动 .C-N 对 称 伸缩 \N-C- 阳 剪 式 振动 、 茶 环 面 内 弯曲 
1487 1487 C-N-C 非 对 称 伸缩 .C- 了 HH 面 外 弯曲 .CH, 面 外 弯曲 CH, 剪 式 振动 
z R 3.2 ARTES e BE E SERS 光谱 分 析 
由 图 3 RI. rH REDE B o FE A TE 983 
i _ 与 1047 cm'' 波 段 处 有 区 别 于 空白 脐橙 样本 的 拉 
2 ož A SREE TITAN EMG E RS HL XT fk 983 与 
= 1041 em- 波段 处 均 有 特征 峰 ， 并 在 1041 em" UE 
段 处 有 6 cm WY, EAT EXE FH TT EAS HEBR 
; 强 的 背景 峰 影响 导 臻 的。 因此， 可 以 将 这 两 处 特 
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tm 图 3 两 组 样本 的 平均 SERS 光谱 
Fig. 3 Average SERS spectra of two groups of samples 


3.3 样本 SERS 强度 检测 条 件 优化 


Fig. 2 Raman Spectrum of theoretical calculation and stan- 


dard substance IMZ 
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的 反应 时 间 不 同 ，SERS 信号 强度 不 同 。 随 着 反 
应 时 间 的 增长 983 和 1047 cm- 处 的 SERS 强度 
均 旦 现 出 持续 下 降 的 趋势 。 产 生 这 种 现象 的 原因 
可 能 是 刚 开 始 吸附 时 ， 脐 橙 果 皮 的 抑 考 哗 与 金 胶 
及 NaBr 溶液 充分 混合 ， 使 得 该 活性 位 点 处 的 
SERS 信 号 得 以 迅速 增强 ""。 吸 附 时 间 过 长 会 
使 金 胶 纳米 粒子 凝聚 过 度 从 而 产生 沉淀 ， 降 低 
抑 霉 唑 与 金 胶结 合 的 效率 “。 但 是 由 于 0 min 
的 反应 时 间 在 实际 操作 中 比较 困难 ， 故 本 研究 
选择 2 min 作 为 最 优 的 反应 时 间 。 


n —1m— 983 cm” 
—e— 1047 cm! 


强度 /au. 


2 
E 


550 


反应 时 间 /min 
图 4 反应 时 间 对 SERS 强度 的 影响 

Fig. 4 Effect of reaction time on SERS intensity 

SERS 主要 是 通过 分 子 吸 附 在 金 和 银 等 粗糙 
金属 的 表面 ， 从 而 使 局 部 表面 产生 等 离 激 元 共 
振 效应 ， 对 SERS 强度 起 到 增强 的 效果 I. gd 
图 5 可 看 出 ， 随 着 金 胶 加 入 量 的 增 大 ，983 与 
1047 cm 波段 处 的 拉 曼 特征 峰 强 度 呈 先 减 小 后 增 
大 再 次 减 小 的 趋势 。 出 现 这 种 现象 可 能 的 原因 是 
当 金 胶 的 量 增 大 到 一 定 程 度 时 ， 有 更 多 的 抑 霉 只 
分 子 吸 附 在 纳米 金 胶 颗粒 的 表面 ， 使 得 983 cm" 
的 特征 峰 在 560 uL HT AY SERS 信号 达到 最 强 , 
1047 cm 的 特征 峰 在 400 uL TRE SERS 信和 号 达到 
最 强 。 而 在 金 腕 加 入 量 为 400 wh WY 983 cm 波段 
处 的 SERS 强 度 为 次 强 点 ， 仅 次 于 560 pL 处 的 
SERS 强度 。 综 合 考 虑 ， 本 研究 选择 最 优 金 胺 加 
人 量 为 400 wh. 

电解 质 溶液 对 SERS 强度 的 增强 效果 取决 于 
离子 对 金 胶 表面 状态 的 影响 及 离子 对 金 胶 纳米 粒 


—m— 983 cm! 
—e— 1047 cnr! 


强度 /a.u. 


iE 


200 上 


Wo — 200 30 4X 50 $9 7 
金 胶 加 入 量 /hL 

图 5 金 胶 加 入 量 对 SERS 强度 的 影响 
Fig. 5 Effects of gold colloid volume on SERS intensity 
子 的 吸附 能 力 所 。 图 6 可 看 出 ， 金 胶 与 CaBr, 的 
混合 溶液 在 983 cm ' 波 段 处 的 SERS 强度 达到 最 
强 ， 金 胶 与 NaBr 的 混合 溶液 在 1047 cm 波段 处 
的 SERS 强 度 最 强 。 出 现 此 现象 的 原因 可 能 是 由 
于 氮 离 子 的 盐 析 能 力 大 于 省 离子 的 盐 析 能 力 ， 导 
致 省 离子 与 金属 粒子 的 表面 络 合 能 力 强 于 握 离 
子 ， 结 果 表 明 含 省 离子 的 化 合 物 (NaBr Ca- 
Br,) XJ SERS 强度 的 增强 效果 明显 大 于 含 氯 离子 
的 化 合 物 (NaCl 和 CaCl,) Xf SERS 强度 的 增强 
效果 中。 其 次 ，NaBr 在 983 cm 特征 峰 处 的 
SERS 强度 仅 次 于 CaBr,。 综 合 考 虑 ， 本 研究 选择 
NaBr 作 为 最 优 的 电解 质 溶液 。 

在 纳米 金 胶 溶液 中 加 入 不 同 量 的 NaBr， 可 
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图 6 不 同 电解 质 溶液 对 SERS 强度 的 影响 
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Fig. 6 Effects of electrolyte solution on SERS intensity 
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能 会 对 SERS 信号 产生 影响 ， 因 此 有 必要 对 电解 
质 溶液 的 加 入 量 进 行 优化 ”。 从 图 7 可 看 出 ， 当 
NaBr 加 入 量 为 10 pL 时，1047 cm 波段 处 的 
SERS 信 号 强度 达到 局 部 最 高 。 当 NaBr 加 入 量 增 
大 到 25 pL 时 ，983 cm PREC AY SERS 强度 达到 
最 高 点 ， 而 1047 cm 波段 处 虽 未 达到 最 强 ， 却 
仅 次 于 最 强 点 。 可 能 的 原因 是 当 讽 族 元 素 的 量 高 
到 足以 聚合 时 ， 使 局 部 的 电磁 场 增强 ， 从 而 增强 
SERS 信和 号。 而 当 加 入 过 量 NaBr 后 ， 会 加 剧 金 纳 
米粒 子 的 聚集 从 而 发 生 沉淀 UU. xti SERS 
言 号 强度 。 因 此 本 试验 确定 了 NaBr 溶 液 的 最 佳 
加 入 量 为 25 uL. 


1200 p 


1000 


400 


200 L1 a 1 n 1 1 " 1 1 : L 
0 5 10 15 20 25 30 


NaBr 加 入 量 /4L 
图 7 NaBr 加 入 量 对 SERS 强度 的 影响 


Fig. 7 Effects of NaBr volume on SERS intensity 


3.4 光谱 预 处 理 


为 减少 光谱 分 析 过 程 中 受到 噪声 、 基 线 漂移 
等 因素 干扰 ， 在 建立 模型 之 前 有 必要 对 原始 光谱 
进行 光谱 预 处 理 ”。 目 前 常用 的 光谱 预 处 理 方 
法 主要 有 air PLS、 归 一 化 、 基 线 校正 、 一 阶 导 
数 、 二 阶 导数 、 标 准 正 态 变 量 (Standard Normal 
Distribution, SNV) 和 多 元 散射 校正 (Multipli- 
cative Scatter Correction, MSC) 等 =, Hp 
线 校正 、 一 阶 导 数 和 二 阶 导数 的 作用 主要 是 扣除 
漂移 或 背景 对 结果 的 影响 ， 在 改变 了 原始 数据 相 
关 性 的 变化 趋势 的 同时 ， 也 有 效 提高 了 某 些 位 置 
光谱 的 相关 性 。 在 光谱 数据 分 析 时 ， 样 本 的 误差 
表现 在 无 规则 分 布 或 样本 颗粒 大 小 各 异 ， 而 


MSC FI SNV 可 以 消除 该 误差 对 光谱 分 析 结 果 的 
影响 ， 还 能 够 提高 原始 光谱 与 预 处 理 后 光谱 之 间 
的 相关 性 。 归 一 化 能 够 消除 由 于 尺寸 差异 过 大 对 
ub usd d 3。 因此 ， 本 研究 以 自 适 
迭代 惩罚 最 小 二 乘法 (Adaptive Iterative Re- 
Penalized Least Squares, air PLS), air 
PLS+ 归 一 化 、air PLS+ 基 线 校 正 、air PLS+ 一 阶 
导数 、air PLS+SNV 和 air PLS+MSC 处 理 后 的 6 
组 光谱 数据 为 研究 对 象 ， 分 别 采 用 这 6 组 光谱 预 
处 理 法 建立 SVR 模 型 并 进行 预测 性 能 的 比较 。 

由 表 2 可 知 ，air PLS, air PLS-SNV, air- 
PLS+MSC 三 组 光谱 预 处 理 方法 所 建立 模型 的 训 
练 集 相 关系 数 (Coefficient of Determination for 
the Training Set, R,) 的 值 均 达到 0.96 以 上 。 且 
三 者 的 训练 集 均 方 根 误 差 (Root-Mean-Square 
Error of Training set, RMSEC) 的 值 相差 不 大 ， 
稳定 在 3.7 附 近 。 然 后 对 该 三 种 方法 的 预测 集 相 
关系 数 (Coefficient ofthe Determinant for the Pre- 
diction Set, R,) 进行 比较 ，air PLS 方法 建立 模 
型 预测 的 预测 集 尺 为 最 大 ， 预 测 集 均 方 根 误差 
(Root-Mean-Square Error of Prediction, RMSEP) 
为 最 小 。 因 此 ， 本 研究 最 终 选 择 air PLS 作为 脐 
TF ze E 4M 2: A SE FA B SVR 预测 模型 的 光谱 预 处 
理 方法 。 

表 2 6 组 光谱 预 处 理 方法 对 SVR 模 型 的 评价 


Table 2 Evaluation of six spectral preprocessing methods 


for SVR model 
预 处 理 R, RMSEC R, RMSEP 
air PLS 0.9620 3.8165 0.9156 4.8407 


air PLS+ 归 一 化 0.7955 9.1730 0.7890 9.0271 


air PLS+ 基 线 校正 0.9105 6.4134 0.9010 5.7491 


air PLS+ 一 阶 导数 0.9461 5.5340 0.8049 8.0074 


air PLS-SNV 0.9702 3.7241 0.6737 9.2194 


air PLS+MSC 0.9751 3.2697 0.6729 9.1564 


3.5 PCA 特 征 提取 


建立 模型 之 前 对 光谱 数据 进行 主 成 分 分 析 
(Principal Component Analysis, n SE 提 
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取 。 一 般 来 说 ， 主 成 分 分 析 时 主 成 分 的 方差 贡献 
率 达 到 9096 以 上 即 可 作为 主要 光谱 信息 的 特 
征 号。 由 表 3 可 知 ， 前 7 个 主 成 分 的 累计 方差 贡 
献 率 已 达 91.94%。 因 此 ， 选 择 前 7 个 主 成 分 的 得 
分 值 作为 SVR 预测 模型 的 输入 值 。 

表 3 前 7 个 主 成 分 因子 的 方差 贡献 率 


Table 3 The variance contribution rates of the first 7 princi- 


pal component factors 


主 成 分 因子 方差 贡献 率 /% 累计 方差 贡献 率 /% 
1 47.52 47.52 
2 17.56 65.08 
3 9.84 74.92 
4 7.52 82.44 
5 4.66 87.10 
6 2.80 89.90 
7 2.04 91.94 


3.6 预测 模型 的 建立 以 及 预测 性 能 分 析 


为 评价 模型 的 预测 精度 和 效果 ， 采 用 了 
SVR、 多 元 线性 回归 (Multiple Linear Regres- 
sion，MLR) 和 偏 最 小 二 乘 回 归 (Partial Least 
Squares Regression，PLSR) 三 种 建 模 方 法 进行 
分 析 比 较 。 其 中 ，SVR 预测 模型 的 类 型 为 epsi- 
lon SVR, {2% KAJ 4 IH BEY eH AL, epsilon (FN 
0.1， 惩 罚 因 子 C 值 为 10。 
通常 选择 R, RMSEP 和 相对 标准 偏差 
(Relative Standard Deviation，RPD) 的 值 来 评价 
模型 性 能 。 其 中 ，R, 和 RMSEP 分 别 用 来 解释 真 
实 值 与 预测 值 的 相关 性 和 精确 度 的 能 力 。R, 越 
大 ，RMSEP 越 小 则 说 明 该 模型 的 拟 合 程度 越 好 ， 
准确 性 越 高 o RPD 值 则 是 对 模型 预测 能 力 的 
评估 RPD 值 越 大 (不 小 于 1)， 说 明 预 测 效果 
越 好 ， 其 预测 结果 越 接近 实测 值 ” 。 

由 表 4 可 得 ， 基 于 SVR 建立 模型 预测 的 尺 ,为 
SVR, MLR, PLSR 三 者 中 的 最 大 ，RMSEP 为 最 
Wy. AIE, SEF SVR SEK XI WE PE Pe TEL I he OR HK 
留 的 预测 值 越 接近 实测 值 ， 能 有 效 提高 模型 预测 
准确 性 。 其 中 ， 最 低 检 出 限 为 1.43 mg/kg (国家 


T (GB 2763-2019 食品 安全 国家 标准 食品 中 农 
药 最 大 残留 限量 》 中 规定 抑 霉 只 在 脐橙 中 残留 最 


大 限量 为 5 mg/kg). 
表 4 不 同 建 模 方法 对 脐橙 果皮 中 抑 老 唑 残留 量 的 预测 
结果 的 影响 


Table 4 Influence of different modeling methods on predic- 


tion results of IMZ residues in navel orange peel 


评价 指标 SVR MLR PLSR 

R, 0.9156 0.9052 0.9064 
RMSEP/( mg:kg ') 4.8407 5.3996 5.2774 
RPD 2.3103 2.3460 1.9830 


为 进一步 做 预测 分 析 ， 在 随机 挑选 35 个 样 
本 建立 模型 之 后 ， 剩 余 的 30 个 样本 作为 预测 集 
进行 验证 ， 即 脐 术 果皮 抑 考 唑 残留 浓度 分 别 为 
38.93、30.45、21.53、14.63、9.24 和 4.05 mg/kg 
的 样本 ， 实 测 的 真实 值 与 预测 值 之 间 的 关系 如 
图 8，R, 为 0.9156，RMSEP 为 4.8407 mg/kg. 


m WAR 
e ”预测 集 
训练 集 拟 合 曲线 


40 | ( 拟 合 方 程 : y=0.8766x+2.3819) 
预测 集 拟 合 曲 线 "s 
( 拟 合 方程 : y=0.8366x+3.4059) "d 


TROU (EL / Cmg: kg) 


0 10 20 30 40 50 
HUSA /mg:kg 
图 8 抑 霉 唑 样本 关系 图 


Fig. 8 IMZ samples relation diagram 
4 结论 


本 研究 探索 了 SERS 技术 在 检测 脐橙 果皮 中 
抑 霉 唑 农药 残留 的 可 行 性 。 首 先 通过 对 比 983 和 
1047 cm"' 两 处 波段 的 SERS gn BE, EFT IM AER IE 
FAHY SERS 检测 条 件 优化 ， 分 别 确定 了 最 优 的 检 
测 条 件 (2 min 的 反应 时 间 ，400 uL 的 金 腕 加 入 
量 ，NaBr 作 为 电解 质 深 液 且 添 加 量 为 253 pL)。 
在 最 优 条 件 下 ， 对 含 抑 霉 唑 的 脐 栖 果 皮 拉 曼 光 谱 
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采用 6 组 不 同 的 光谱 预 处 理 方法 ， 分 别 对 这 6 组 
光谱 所 建立 模型 的 预测 性 能 进行 对 比 ， 最 终 采用 
air PLS 的 光谱 预 处 理 方法 。 经 过 PCA 的 特征 提 
取 ， 选 择 前 7 个 主 成 分 得 分 作为 脐橙 果皮 抑 霉 唑 
残留 的 SVR 预测 模型 的 输入 值 。 最 后 SVR 预测 
模型 的 R, 可 达 0.9156，RMSEP 为 4.8407 mg/kg, 
RPD 为 2.3103。 由 于 本 研究 所 采用 的 计算 实际 残 
留 浓度 的 方法 为 称 重 法 ， 挥 发 时 间 不 一 致 容易 造 
成 称 量 误差 ， 进 而 影响 最 终 实际 残留 浓度 的 计 
算 。 因 此 实际 残留 浓度 的 计算 方法 还 有 待 进一步 
改善 。 研 究 结 果 表 明 ，SERS 技 术 可 实现 脐橙 果 
皮 抑 霉 唑 残留 的 快速 检测 。 
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Rapid Detection of Imazalil Residues in Navel Orange 
Peel Using Surface-Enhanced Raman Spectroscopy 


ZHANG Sha", LIU Muhua", CHEN Jinyin’, ZHAO Jinhui 一 


(1. College of Engineering, Jiangxi Agricultural University/Key Laboratory of Modern Agricultural Equipment in Jiangxi 
Province, Nanchang 330045, China; 2. Jiangxi Provincial Collaborative Innovation Center of Key Technologies and Quality 
and Safety in Post-Harvest Processing of Fruits and Vegetables, Nanchang 330045, China) 


Abstract: Imazalil, a preservative for navel orange in the process of postharvest processing, is easy to seep into the flesh through 
the peel and produce residues in the flesh, which is vulnerable to cause endanger to human body if it was eaten accidentally. 
Base on this, a fast detection method of imazalil residues in navel orange peel ,namely surface-enhanced Raman spectroscopy 
(SERS) was proposed in this study. Firstly, the SERS detection conditions of imazalil residues in navel orange peel were opti- 
mized, and the optimal detection conditions were determined as follows: Reaction time of 2 min, gold colloid of 400 uL, NaBr 
as electrolyte solution, NaBr dosage of 25 uL. Based on the above optimal conditions, 6 groups of spectral data processed by 
adaptive iterative penalized least squares (air PLS), air PLS combination with normalization, air PLS combination with baseline 
correction, air PLS combination with first derivative, air PLS combination with standard normal distribution (SNV), air PLS 
combination with multiplicative scatter correction (MSC) were used to establish support vector regression (SVR) models and 
compare the models prediction performance. And air PLS method was selected as the spectral pretreatment method, because the 
value of correlation coefficient computed value of prediction set (R,) is the largest, and the value of root mean square error cal- 
culated value of the prediction set (RMSEP) is the smallest. Then, principal component analysis (PCA) was used to extract the 
features from spectral data, and the first seven principal component scores were selected as the input values of SVR prediction 
model. SVR, multiple linear regression (MLR) and partial least squares regression (PLSR) were used to analyze and compare 
the prediction performances. The R, value of prediction set of SVR prediction model could reach 0.9156, the RMSEP value of 
their prediction set was 4.8407 mg/kg, and the relative standard deviation computation value (RPD) was 2.3103, which indicat- 
ed that the closer the predicted value of imazalil residue on navel orange surface based on SVR algorithm was to the measured 
value, the more effective the prediction accuracy of the model could be. The above data indicated that the speedy detection of 
imazalil residues in navel orange peel could be emploied by SERS coupled with PCA and SVR modeling method. 

Key words: navel orange; imazalil; surface-enhanced Raman spectroscopy; support vector regression; multiple linear regression; 


partial least squares regression 
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